Superconductivity in mesoscopic high-T c superconducting particles 
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Based on the Hubbard model in the framework of non-phonon kinematical mechanism and 
taking into account the discreetness of an electronic energy spectrum, the superconducting critical 
temperature of a mesoscopic high-T c sphere is analyzed as a function of doping and as a function of 
a particle's radius. The critical temperature T c is found to be an oscillating function of the radius 
of a particle. The size-dependent doping regime is revealed in high-T c nanoparticles. Our analysis 
| shows that each oscillation in T c corresponds to the increase of a number of the energy levels in the 

sphere by one. The amplitude of oscillations of T c increases with decreasing R and can reach a value 
of 6 K for nanopartilces with sizes about 25 nm, in a good agreement with experimental studies of 
YBa2Cu307_,5 nanoparticles. 
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£3 ' I. INTRODUCTION 

After experiments by Ralph et al. ' on the electronic states in nanosize Al particles there has been increasing 
experimental and theoretical interest in properties of ultrasmall superconducting samples with discrete fermionic 
energy spectrum (see 3,4 and Refs therein). Incorporation of the nanosize particles as a part of single-electron tunneling 
transistor schemes allowed for measurements of excitation spectra as a function of the number of electrons in a 
superconducting particle. The large spectroscopic gap between the lowest energy level and all the others which can 
be driven to zero by application of magnetic fields has been explained by superconducting pairing. 

Till now experiments have been completed for the low T c elemental superconductors (Al, Sn, ...), showing traces 
of superconducting pairing for metallic particles even with the electron energy spacing being larger than the super- 
conducting energy gap: d > A. For elemental superconductors in the order of magnitude A ~ T c , the energy spacing 
is d ~ h /(m e R 2 ) and for these nanoparticles their size R < h/y/m e T c . Thus, to satisfy the regime d > A for charac- 
teristic critical temperature values T c ~1K the nanoparticle size should not exceed ~ 30 nm 5 . As shown by Tinkham 
team, the problems of the sample preparation with these parameters have been solved by modern technology. In such 
y—\ . nanoparticles from elemental superconductors their average size is smaller or comparable with the coherence length 
t— I ' and fluctuations of superconducting order parameter do not allow the mean- field treatment. On this way the recent 
theoretical work includes studies of mesoscopic systems with equally spaced levels 6 , Wigner-Dyson level spacing 7 , 
parity effect 8 , fluctuation effects 9 , fixed-N canonical treatment 10 and the suppressing role of the non-pair portion of 
the interaction 11 . The differences in energy level spectra for even and odd numbers of electrons are being interpreted 
as a result of superconducting paring interactions. 

But today the nanotechnology is on the threshold of production of nanosize samples from high-T c cuprates 12-14 . 
More recently, a several degrees enhancement of the onset transition temperature T c has been observed in super- 
conducting dots formed by the high-Tc superconductor YBa 2 Cu307_,5 (YBCO) 15 and in YBa2Cu307_5 films with 
detected spherical particles on their surface 15 . The dots were formed when powder YBCO particles with a diam- 
O ■ eter < 600 nm were exposed to a RF plasma to produce a Coulomb crystal. The plasma damaged the particles 
J> \ and caused 20-25 nm sized isolated islands of the correct stoichiometry to segregate within each particle. Possible 
technological applications of superconducting nanoparticles make it desirable, to engineer their energy spectra and 
the pairing characteristics. The crucial question that has arisen in this connection is: how can the superconducting 
critical temperature for a high-Tc nanoparticle be controlled by changing its size and the doping degree? 

In nan opar ticles made from high-T c materials it's questionable to reach the regime d > A corresponding to size 
R < h/VmA = h\/2/^mT c (2A/T c ). Really for high-T c superconductors the BCS ratios 2A/T C are quite large (9 
for La2-a;Sra;Cu04, 6-8 for YBa2Cu307_5, and even 12 for Bi2Sr2- a; La a ;Cu06 - see 16 and Refs. therein) and the 
magnitudes of an effective electron mass m are enhanced due to narrow energy bands of high-T c materials. As a 
result, the high-T c nanoparticle size R is becoming comparable with a unit cell size (it's lattice constant parameter 
c ~ 1.2 nm in YBCO) or even less than the latter with a significant deviations of a chemical composition and the 
intra-size vibrations from those in a bulk mother substance. Physically, it means that for the nanosize particles 
formed from stochiometric high-T c material we can not reach regime A < d. Thus, the field of the nano-size high-T c 
superconductors is renewing the earlier Anderson idea 17 about the disappearance of superconductivity in particles, 
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which are so small that the "granularity" of an electron energy d is larger than the superconducting energy gap A. 
So, for the nanosized high-T c particles we arc returning back to the limit d < A. 

Long time ago the studies of the superconducting-transion temperature T c of metallic films of Al, Sn, In, Zn 
prepared by vacuum deposition onto substrate held at cryogenic temperatures revealed the quantum size effect of 
the T c -oscillations with a film thickness 18-24 (see also review 25 ). In the simplest approaches 26-28 the quantum size 
effects follow from oscillations in the electronic density of states with thickness which in turn lead to oscillations with 
thickness of T c and other film quantities. 

Here in the Anderson regime d<Awe consider a similar mechanism of finite-size effects aiming the goal to apply 
it to realistic high-T c nano-particles. Worthy to note, that an evaluated size of nanoparticles under consideration (see 
next section) exceeds the small coherence length scale with neglect of fluctuations and the Aslamazov-Larkin effect 29 . 
The method is based on the study of the non-phonon kinematical superconductivity (see 30 and references therein) 
in systems with strongly interacting electrons and an intrinsic "granularity" of energy, induced by the discreteness 
of electron energy levels in a nanoparticlc. Formally the non-phonon mechanism allows to avoid difficulties caused 
by phonon cut-off energies in small particles. In the model we assume the boundary condition on the electronic 
wave- functions is that they vanish on the surface of spherical nanoparticlc and we neglect the surface effects. 



II. THE MODEL 



In the present communication we tackle the problem starting from the strongly interacting limit of the Hubbard 
model (U » t), H = — ^ZtijXf X®' 7 — [iJ2 n ii f° r correlated electrons (fi: their chemical potential, nf. electron 
density operator on site, tif intersite electron hopping) in the YBCO nanoparticles of a spherical shape. As usual, 
X 0rT denote a projection Hubbard operator with an index "Oct" labeling an intrasite transition from an energy level 
occupied by an electron with spin projection ct to an empty one: a — > 0. In high-T c materials the characteristic strong 
electron interactions remove the orbital degeneracy of energy levels. The electron energies have been extracted 30 ' 31 
from zeros of the inverse Green's function, obeying the Dyson equation 

V- 1 (iu) = [V°{iuj)]- 1 +Uj, (1) 

where the first order self-energy is itself an electron hopping t t j inside a spherical nanoparticlc The zeroth order 
Green's function is defined as T>°(t) — — (TX 0(7 (T)X a0 (0)} or in Matsubara u) = (2n+ l)nT -representation it is such 
as 

V°(iuj) = - . (2) 

— KJ n - /I 

The correlation factor, 

/ = (A 00 + X aa ) = 1 - (X™) (3) 

n 1 + x 
= l-n w =l-- = —, 

is governed by an average electron density per site n or/and by a doping parameter x = 1 — n (0 < x < 1). 

In the momentum representation of bulk material the correlated electron energy is £ = ft p — fi with an energy 
dispersion values t p located inside the band: — w < t p < w. The superconducting instability is driven by the 
homogeneous Bcthe-Salpeter equation for a vertex T of two scattering fermions with opposite spin orientations and 
opposite momenta in their reference system: 

r = 2Tj2tpG-A-p)GApW(p), (4) 

where G ul {p) = (—iuj + ^p)^ 1 is the normal Green's function. After summation over Matsubara frequencies this 
equation can be converted to equation 
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represented in momentum representation. 
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In the experiment 15 both the initial YBCO material and nanosized (20-25 nm) islands produced by plasma damage 
are polycrystalline particles with chaotic distribution of crystallographic directions. Therefore, we can model each 
nanosized island by isotropic sphere of radius R. As in the case of dirty superconductors 17 , in a spherical particle 
of radius R with the quantized electronic energy levels, E v = h 2 v 2 / (2mi? 2 ), the momenta are not good quantum 
numbers. In the absence of non-elastic processes the suitable quantum numbers are the eigen-energies of the problem, 
shifted by a half-bandwidth w: e v = E v — w (—w < e„ < w). For HTSC's of interest such as YBa 2 Cu307_,5 a 
half-bandwidth parameter has characteristic magnitude w ~ 0.5 eV 32 . The lowest energy level is at the bottom 
of the starting bulk energy band, e = — w, whereas the highest one, e Vo , coincides with the top of the band: 
= [Ry/2mw/h\ ([...] stands for an integer number). So, in an energy representation the superconducting critical 
temperature T c is determined by the following equation: 

•, "o El-\ / (Ex. -I) f - Ji.\ 

u—l V w ) •> w \ w / 

where the summation is running over quantized energy levels till the upper level ^ ,' the chemical potential \i has a 
meaning of the highest occupied energy level. The doping parameter x is defined by 

l-x = 2T J2 e iuJn5 V{u, uj n ). (7) 

TIM — 1 

From here one can find the equation for the chemical potential ft: 

(8) 



1 + x 




The system of equations, Eqs. (6), (8) should be solved self-consistently to obtain T c . 

The applied arithmetics is valid until the superconducting pairing is allowed for electrons belonging to different 
energy levels in a nanoparticlc, i.e. the thermal energy should exceed the distance between the discrete energy levels, 

namely T c >h 2 v 2 — (i/q — l) 2 /(2mi? 2 ), wherefrom it follows the restriction for the high-T c nanoparticle radii in our 

approach: 

R > Rq = — — \ / — . 

For estimate we assume the effective mass enhancement m/m e w 10 from the earlier optical reflectivity measurements 33 
for the YBCO bulk crystals. Noteworthy, the more recent data for YBCO 34 give closed estimate m « 8 — 9m e under 
the assumption of a twice band mass enhancement due to electron correlations: = 2m e . Such estimate m w 10m e 
is also in an agreement with theoretical evaluations 35 ' 36 for an effective electron mass in YBCO. Finally, for the YBCO 
nanoparticle of interest 15 with T c — 92 K and a half-bandwidth w = 0.725 cV (this value of w corresponds in our 
calculations to the "bulk" critical temperature T c — 92 K) one can conclude the lower critical radius of a spherical 
YBCO nanoparticle such as Ro = 15 nm to apply the theory. The nanoparticle with a critical radius still contains 
757T • 10 3 unit cells of YBCO (the YBCO unit cell volume is <~ 0.18 nm 3 ), which is in a reasonable agreement with a 
neglect of surface and fluctuation effects. As it will be seen from numerical results in next sections, although we are 
working in the bulk superconducting regime, the quantum size effects due to a " granularity" of an energy in high-T c 
nanosize particles are still traced. 



III. RESULTS AND DISCUSSION 



A. The critical temperature T c of a mesoscopic high-T c superconducting sphere as a function of doping x 

Based on Eqs. (6) and (8), we analyze the critical temperature of a mesoscopic high-T c superconducting sphere as 
a function of doping. We have performed calculations for various radii R of a superconducting sphere varying from 
R = 15 nm to 1000 nm. In Fig. 1, the critical temperature T c of a mesoscopic high-T c superconducting sphere is 
plotted as a function of doping x for particles with radii R — 15 nm, 16 nm, 20 nm, 30 nm and 1000 nm. For any 
radius within the above region, T c as a function of x is characterized by the following behavior: T c is zero until some 
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minimal doping x m i n , which provides a particle to be superconducting; then T c increases (underdoped regime) and 
reaches a maximum at some value of doping (optimal doping) ; further increase of x leads to decrease of T c (overdoped 
regime) (see Fig. 1). 
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FIG. 1. The critical temperature T c as a function of doping x, for a mesoscopic high-T c superconducting spherical particle. 
The radii of the particles are R = 15 nm (a), 16 nm (b), 20 nm (c), 30 nm (d) and 1000 nm (e). The value of the bandwidth 
parameter for YBa2Cu307-<s is w = 0.725 eV. The inset shows upper parts of the curves T c = T c (x) for temperatures ranging 
in the interval from 60 to 100 K. 

The maximal value of T c — T c (x) is a non-monotonous function of R. For example, at R = 15 nm the maximum 
of T c is 86.6 K, whereas at R = 16 nm the maximumal value of T c is as high as 96.1 K. With increasing radius, a 
deviation in the maximal value of T c reduces, and for large radii (more than 100 nm) the maximal values of T c are 
close to each other and lie in the vicinity to T c w 92 K. A further increase of the radius of the sphere does not lead 
to any appreciable changes of T c . 

A position of the maximum of T c = T c (x) itself also varies as a function of radius R. For R — 15 nm, the maximum 
of T c — T c (x) is reached at a higher x as compared to "bulk" value (R = 1000 nm). On the contrary, for R = 16 nm, 
the maximum of T c = T c (x) is at a lower x as compared to the "bulk" value. 

The critical temperature of a mesoscopic high-T c superconducting sphere as a function of doping is characterized 
by some minimal value of doping and by oscillating behaviour in strongly overdoped regime (see Fig. 1). Although 
the critical temperature of a mesoscopic high-T c superconducting sphere is itself a smooth function of doping, the 
minimal doping x m i n occurs to be an oscillating function of the radius of the particle. The minimal doping, which 
provides the particle to be superconducting, is very sensitive to the bandwidth parameter w. Varying the parameter 
w, we can conclude that reduction of the bandwidth leads to increase of the average value of x m i n and to decrease of 
the period of oscillations. 



In this section, we analyze in detail the oscillations in T c that were noted in the previous subsection. The critical 
temperature T c is calculated as a function of the radius R of the sphere for various values of doping x. In Fig. 2, 
T c = T C (R) is shown for x = 0.2 (that is close to the optimal doping). The critical temperature T c as a function of 
radius R is characterized by an oscillating behavior. The amplitude of these oscillations decreases with increasing 
R. It is maximal at small R: T c = T C (R) varies from about 86.4 K up to 97.8 K in the interval from R = 15 nm to 
18 nm (see the upper inset to Fig. 2). Then, as shown in the bottom inset to Fig. 2, the amplitude of the oscillations 
decreases very fast up to R « 100 nm. For larger R, T c = T C (R) is characterized by a slow decrease versus R. The 
function T c — T C (R) "saturates" with increasing R and reaches the "bulk" value 92 K. 



B. The critical temperature of a mesoscopic high-T c superconducting sphere versus its radius 
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FIG. 2. The critical temperature T c as a function of the radius J? of a mesoscopic high-T c superconducting particle, at the 
optimal doping x = 0.2, for radii from R = 15 nm to 1000 nm. In the bottom inset: T c as a function R for radii from R = 15 nm 
to 75 nm. The top inset: T c = T C (R) for R — 15 nm to 18 nm. Each oscillation in T c = T C (R) corresponds to the change of 
the number of levels in the system by one. 

The oscillations in the function T c — T C (R) are due to the discrete density of states in the superconducting sphere. 
Recall that we analyze high-T c spheres with a confinement of carriers in all three dimensions. Our analysis shows that 
each oscillation in T c with increasing R corresponds to the increase of a number of the energy levels in the sphere by 
one. 

In the overdoped regime, the function T c = T C (R) is analyzed for x = 0.225. The results of the calculations are 
shown in Fig. 3 for R = 15 nm to 1000 nm The "bulk" value of T c for x = 0.225 {T cMlk = 62 A K) is much lower than 
that for the optimal doping. Therefore, slowly varying the doping parameter x, we can model various real samples 
used in experiments. 
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FIG. 3. The critical temperature T c as a function of R in the overdoped regime, x = 0.225. The range of particles radii and 
the bandwidth parameter are the same as in Fig. 2. 

In the paper by Menon et al. 15 , a slight enhancement in the critical temperature T c was detected for YBaaCu307_i5 
plasma-damaged nanoparticles with sizes about 25 nm. To find the transition temperature in the particles, their dc 



5 



susceptibility was measured as a function of temperature in a SQUID magnetometer. The transition temperature for 
the onset of superconductivity is defined as the temperature at which the susceptibility begins to drop 15 . In some 
cases, the increase in critical temperature was only 1 K (Tc.btiifc = 92 K), whereas the increase reached 6 K for some 
other samples 15 . From our calculations it follows that, for the optimal doping (characterized by the "bulk" value 
Tc.buik = 92 K) the enhancement of the critical temperature for nanoparticles with the radii R = 15 nm to R = 18 nm 
can be up to 6 K (see upper inset to Fig. 2). This is in a rather good agreement with the experimental results 15 . 
Moreover, the revealed oscillations of T c as a function of R explain the non-regular character of the enhancement of 
T c in different samples 15 . 

IV. CONCLUSIONS 

The superconducting critical temperature of a mesoscopic high-T c sphere has been analyzed as a function of doping 
and as a function of a particle's radius. The approach is based on the Hubbard model as the most fundamental model 
of the electronic interactions in high-T c materials. Quantum Monte Carlo studies of this model (see 37 and references 
therein) concluded that the enhanced superconducting correlations are not predominant at low temperatures or in the 
ground state. Although some of numerical works supported the existence of superconducting pairing in the Hubbard 
model 38,39 , these quantum Monte Carlo computations have restricted parameter space such as the system size, the 
doping, the low values of the Hubbard energy etc. The modern Monte Carlo 40-42 and other 43 metods allow one to treat 
a wider parameter space with large magnitude of the Hubbard energy and reveal the superconducting correlations in 
the ground state of the two-dimensional Hubbard model. So, the numerical simulation of superconductivity in the 
ground state of the Hubbard model still remains the issue for the serious problem of the two-dimensional nature of 
cupratc high-T c superconductors. As distinct from that, the kinematical mechanism (see Ref. 30 and references therein) 
considers the superconducting instability of strongly correlated electrons in normal state beyond the knowledge of their 
ground state properties. Therefore the present study of superconductivity is based on the non-phonon kinematical 
mechanism in the framework of the Hubbard model taking into account the "granularity" of electron energies (the 
discreteness of an electronic energy spectrum) in nanoparticles. The chaotic crystallite structure allowed us to model 
nanoparticles by the isotropic spherical particles. This point made suitable the applied three-dimensional scenario of 
the kinematical attraction for YBCO nanoparticles under study. Worthy to note that in one-band two-dimensional 
Hubbard model the kinematical superconducting instability occurs at bulk doping x < 1/3 as in variational Monte 
Carlo method (see, e.g. , Ref. 40 ) . Stimulated by the recent successful results on technology of a nano-samplc preparation 
from YBCO 12-15 , the typical superconducting compound YBa2Cu307_5 has been taken as a mother substance forming 
the nanoparticle. We have shown that the critical temperature of a mesoscopic high-T c superconducting sphere is 
characterized by a non-monotonic dependence on doping as in bulk material. However, contrary to the bulk high-T c 
materials the T c -maximum (corresponding to the optimal doping) is shifted to the lower doping regime. A bulk 
material with a decreased T c value in overdoped regime can form a nanoparticle with a maximal T c value which 
means that it is in an effective optimally doped regime. Also an optimally doped bulk high-T c material behaves as an 
undcrdoped substance in nanoparticles. The revealed size-dependent doping regime in high-T c nanoparticles is new 
effect in comparison with quantum size effects in elemental low-T c nanoparticles. 

Although we worked in the Anderson regime d < A, i.e. neglecting the parity effects, and due to the small 
coherence length of high- T' c s, £ < i?o = 15 nm, neglecting the fluctuation effects, it have been established that the 
critical temperature of a mesoscopic high-T c superconducting sphere is an oscillating function of a large sphere size 
R > i?o- For sphere radii with value R > 1000 nm and doping x = 0.2 — 0.25 the T c -oscillations actually disappear. 
The oscillating quantum size effects appear whenever a size-dependent energy level passes through the chemical 
potential as the size of a nanosphcrc is varied. At each resonance a new quantized energy level starts to contribute 
to an electron density of states and superconducting quantities. For chemical potential nearly zero (x <~ 1/3) this 
oscillating behaviour is more pronounced and results in new T c oscillations with doping in strongly overdoped regime 
(cf. Fig. 1). The amplitude of T c -oscillations is enhanced for lower radii of mesoscopic spheres. For an electron 
confined to a sphere of radius R, the radial part of the wave function is proportional to j\{kR), i.e. to sm(kR)/kR 
for electrons without angular momenta. Therefore the boundary condition yields the relationship kR = nir with a 
period of oscillations such as AR = ir/k = Xf/2, a half of de Broglie wave-length. The calculations are performed 
for spherical YBCO nanoparticles which radius exceeds the critical one i?o = 15 nm. The nanoparticle with a critical 
radius Ro still contains 75ir ■ 10 3 unit cells of YBCO (the YBCO unit cell volume is <~ 0.180 nm 3 ). For the optimally 
doped regime, the enhancement of the superconducting critical temperature for the particles with sizes R = 25 nm 
can be up to 6 K what is in agreement with the experimental findings for YBa2Cu307_,5 nanoparticles in Ref. 15 . 
Also we can not exclude the visualized enhancement AT C >1K in YBa2Cu307_,5 films in Ref. 13 due to detected 
spherical particle formation on their surface in the process of laser ablation with increased deposition temperatures 
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> 790° C. The developed theoretical approach can be applied to nanoparticles prepared from other high-T c materials 
with narrow electronic energy bands in bulk. 
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